Defect solitons supported by nonlocal PT symmetric superlattices 
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The existence and stability of defect sohtons supported by parity-time (PT) symmetric super- 
lattices with nonlocal nonlinearity are investigated. In the semi-infinite gap, in-phase solitons are 
found to exist stably for positive or zero defects, but can not exist in the presence of negative defects 
with strong nonlocality. In the first gap, out-of-phase solitons are stable for positive or zero defects, 
whereas in-phase solitons are stable for negative defects. The dependence of soliton stabilities on 
modulation depth of the PT potentials is studied. It is interesting that solitons can exist stably for 
positive and zero defects when the PT potentials are above the phase transition points. 
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I. INTRODUCTION 

Over the last few years, people have paid much at- 
tention to parity-time (PT) symmetric systems 
Bender et al have shown that a family of Hamiltonian 
with PT symmetry can still exhibit entirely real spec- 
tra, despite the non-Hermitian Hamiltonian [T(^ . PT 
symmetry requires that the complex potential satisfies 
V(x) = V* {—x). For a PT symmetric Hamiltonian, there 
exists a critical threshold above which its eigenvalues are 
not real but become complex, and the system undergoes 
a phase transition because of the spontaneous PT sym- 
metry breaking flo'|. Quite recently the notion of PT 
symmetry was introduced within the framework of op- 
tics This suggestion was based on judicious designs 
that involve both optical gain /loss regions and the pro- 
cess of index guiding pi [sl Unusual properties such 
as unidirectional invisibility [l2| and nonlinear switch- 
ing structures (l3l [T^ have been found in PT symmetric 
structures. PT symmetry breaking in synthetic optical 
systems has been observed experimentally in semiconduc- 
tor heterostructure Q, photorefractive crystals [6], and 
ii?C circuits 0- 

Soliton phenomena in optical lattices and superlattices 
have attracted considerable attention [T5l - [2l| . Optical su- 
perlattices can be considered as regular lattices hosting 
a periodic chain of defects. Superlattices are fascinating 
because the structures exhibit collective properties not 
shared by either constituent, and these properties can 
be controlled through variation of the structural param- 
eters [l^l. Defect solitons in PT symmetric local lattices 
and superlattices have been studied and stable solitons 
are not found in the first gap for positive or zero defects 
[l^ll^. However, defect solitons in PT symmetric super- 
lattices with nonlocal nonlinearity have not been studied. 
It is noteworthy that the nonlinearity in photorefractive 
media, in which Riiter et al have observed the PT sym- 
metric wave propagation [6] , is nonlocal nonlinearity |25l - 
[27} . Nonlocality pla ys a n important role in many areas 
of nonlinear physics [28| , drastically modifies the proper- 
ties of solitons and improves the soliton stability [28l - l30| . 



Therefore it is worthy to study the properties of solitons 
in PT symmetric optical superlattices with nonlocal non- 
linearity. 

In this paper, we study defect solitons supported by the 
PT symmetric superlattices in nonlocal nonlinear media. 
We find that the nonlocality expands stability ranges of 
solitons, especially in the first gap. Both out-of-phase 
and in-phase solitons can be stable in the first gap for 
positive, zero, and negative defects. In-phase solitons can 
exist stably in the semi-infinite gap for positive or zero 
defects, but exist unstably for negative defects. When the 
modulation depth of the PT potentials is small, solitons 
can exist stably for positive and zero defects, even if the 
PT potentials are above the phase transition points. 



II. MODEL 

We consider the propagation of light beam in PT sym- 
metric defective superlattices with Kerr-type nonlocal 
nonlinearity. The evolution of complex amplitude U of 
the light fields can be described by following dimension- 
less nonlinear Schrodinger equation. 



.dU d^U 



+ + [Vix) +iWix)]U + nU = 0, (1) 



dz dx^ 



dx^ 



(2) 



where x and z are the transverse and longitudinal coor- 
dinates, respectively, n is the nonlinear refractive-index 
change, d stands for the degree of nonlocality of the non- 
linear response. This type of nonlinear response with a 
finite region of nonlocality exists in many real physical 
systems, for instance, the photorefractive crystals used in 
the experimental observation of PT symmetry in optics 
0. When d — >■ 0, Eq. ^ reduces to n = and above 
equations reduce to the local case. V{x) and W{x) are 
the real and imaginary parts of the PT symmetric defec- 
tive superlattices, respectively, which are assumed in this 
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paper as 



V{x) = Vb[ei cos2(x) + (1 - ei) sm^(2x)] 

X [l + eexp(-a;**/128)], 
W{x) = Wo8in{2x). 



(3) 
(4) 



Here e and ci represent the strength of the defect and 
the modulation parameter of superlattice, respectively. 
When 0.1 < ei < 0.7, Eq. ((21) has the superlattice 
shape fis'l. Without losing of generality, we take ei = 0.5 
throughout the paper. The parameters Vq and Wq rep- 
resent the modulation depth of the real and imaginary 
parts of the PT symmetric potentials, respectively. 
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FIG. 1. (Color online) (a) Band structure of the superlattice 
with Vb = 6, Wo = 0.4 and ei = 0.5. (b) Band structure 
of the superlattice with different Wo values, where Vo = 6 
and ei — 0.5. (c) The phase transition point Woe versus 
Vo for different ei values, (d)-(f) Lattice intensity profiles 
for the superlattice with defects (d) e — 0.3, (e) e = 0, and 
(f)e — —0.3 Blue solid lines and red dotted lines represent the 
real and imaginary parts, respectively. 

The PT symmetric superlattices in Eq. ([T]) have a 
Bloch band structure when e = and n = Q. The Bloch 
band structure obtained by the plane wave expansion 
method for ei = 0.5, Vb = 6 and Wq = 0.4 is shown 
in Figllja). One can see that the region of the semi- 
infinite gap is 6 > 3.27, and the first and second gaps 
are 1.56 < b < 2.77 and -1.88 < b < -0.19, respec- 
tively, where 6 is the propagation constant. Figure [IJb) 
shows the Bloch band structure for Vb = 6 with different 
Wq values. The region of the first gap decreases with 
increasing Wq. The first gap disappears at Wq — 1.035, 
which is the phase transition point Wqc- When the sys- 
tem is above the phase transition point, i.e. Wq > 1.035, 
the band structure becomes complex. We find that, for 
the PT symmetric superlattices, Wqc changes with Vq 
nonlinearly, as shown in Fig. [TJc). However, for the 
ordinary lattices(ei = l), the Wqc increases linearly with 
increasing Vq, i.e, Wqc/Vq =0.5 [H, Ql- In this paper, we 
study the PT symmetric superlattices below or above the 
phase transition points, respectively. Figures [Hd)-[IJf) 



show the intensity distributions of the PT superlattice 
potentials with the strength of defects e = 0, e = 0.3 and 
e = —0.3, respectively, e = corresponds to the uniform 
superlattice. 

We search for stationary solutions to Eqs. ([1]) and ([2]) 
in the form U = /(x) exp(i6z), where f{x) is a complex 
function and satisfies equations. 
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The solutions of defect solitons are gotten numerically 
from Eqs. (0 and ^ by the modified squared-operator 
method 31] and shown in the next section. To eluci- 
date the stability of defect solitons, we search for the 
perturbed solutions to Eqs. ((Ij and ([2]) in the form 
U{x,z) = [f{x) + u{x,z) + iv{x, z)]e'}q){ibz), where the 
real [u{x, z)] and imaginary [v{x, z)] parts of the pertur- 
bation can grow with a complex rate S upon propagation, 
i.e. u(x,z) = u{x) exTp{i6z) and v{x,z) = v{x) expiiSz). 
Substituting the perturbed solution into Eq. ([1} and lin- 
earization of it around the stationary solution f{x) yields 
the eigenvalue problem 



Sv 



Su ■ 



+ (n - b)u + (Vu - Wv) 
+ Reifix)] JZ, 2G(x - OniOMfiOn 
+ Re[f{x)] 2G{x - OlmifiOHOd^, (7) 
- - (n - 6)w - {Wu + Vv) 

- Im[f{x)] J^^ 2G{x - e}Im[fmv{m 

- Im[J{x)] 2G{x - OReifiOHOdt (8) 



Here G{x) = [l/{2Vd)]exp{-\x\/Vd) is the response 
function of the nonlocal media. Above eigenvalue prob- 
lem is solved numerically by the original-operator iter- 
ation methodfl^l to find the maximum value of Re{6). 
If Re{S) > 0, solitons are unstable. Otherwise, they are 
stable. 



III. NUMERICAL RESULTS 

In the nonlocal PT symmetric defective superlattices, 
we find two types of defect solitons. The first type is 
in-phase soliton, whose real part of optical fields is sym- 
metric and the imaginary part is antisymmetric. The 
other type is out-of-phase soliton, whose real part is an- 
tisymmetric and the imaginary part is symmetric. It is 
noteworthy that if U — f{x) exp(ibz) is a solution to 
Eqs. [1] and [U we also have a series of solution, i.e. 
U — f{x)e^^exp{ibz), where 9 is an arbitrary initial 
phase. We know that f{x) and f{x)e^^ represent the 
same physical soliton despite the initial phase, but their 
parities are uncertain. Due to PT symmetry, one can 
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FIG. 2. (Color online) (a) Soliton power P versus propagation 
constant 6; (b) Cutoff point be. versus the depth of defects, (c) 
Unstable growth rate Re[5) versus propagation constant h for 
out-of-phase solitons in the first gap. For all cases Vb = 6, 
Wo — 0.4, and e — 0.3 except for (b). 
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FIG. 3. (Color online) The complex fields (solid blue: 
real part; dotted red: imaginary part) and refractive-index 
changes (dashed black) for in-phase solitons in the semi- 
infinite gap at (a) 6 = 5, d = 1, and (b) fe = 5, d = 3, 
and for out-of-phase solitons in the first gap at (c) b = 2.2, 
d = 1, (d) 6 = 2.2, d = 3, (e) 6 = 2.6, d = 1, (f) & = 2.6, d = 3, 
respectively. For all cases e = 0.3, Vb = 6 and Wo = 0.4. 



find a proper initial phase to guarantee that both the real 
and imaginary parts of solitons have fixed parity (odd or 
even). When we use a real function as the initial trial 
in the iteration program solving Eqs. ([5]) and ([6]), the fi- 
nal convergency solution has a certain parity and its real 
part is always bigger than its imaginary part. Then we 
can distinguish the in-phase and out-of-phase solitons in 
this occasion. We find that in the first gap, out-of-phase 
solitons can exist stably for positive and zero defects, 
whereas in-phase solitons can exist stably for negative 
defects. In the semi-infinite gap, we find in-phase soli- 
tons for both positive, zero, and negative defects. 

For positive defects, we assume e — 0.3 and the re- 
sults are shown in Figs. [^JHl We find in-phase solitons 
in the semi-infinite gap and out-of-phase solitons in the 
first gap, respectively. Figure [2ja) shows that the power 
of solitons [defined as P = \f{x)\'^dx] for both in- 
phase and out-of-phase solitons increases almost linearly 
as increasing of propagation constant. There exists a cut- 
off point of propagation constant, below which the defect 
solitons vanish. The propagation constants of the cutoff 
points do not depend on the nonlocality degree. This fea- 
ture is similar to the case of traditional uniform lattices 



FIG. 4. (Color online) (a)-(f) Evolutions of defect solitons 
corresponding to Figs. E^a)- EJ^f), respectively. 



in nonlocal media [291 . The propagation constants of the 
cutoff points for both in-phase and out-of-phase solitons 
increase with increasing the depth of positive defects, as 
shown in Fig. ^h). When e > 0.725, the out-of-phase 
solitons can not exist in the first gap. 

We find that in-phase solitons are stable in the whole 
regime where solitons exist in the semi-infinite gap when 
Wo = 0.4, like its counterparts in local FT lattices. Out- 
of-phase solitons are found to exist stably for the low 
propagation constants in the first gap, as shown in Fig. 
I^c) . The stability regions of defect solitons decrease with 
increasing Wq. In-phase solitons in the semi- infinite gap 
become unstable for large propagation when Wq > 0.46 
and their stability region vanishes totally as Wq > 0.58. 
Out-of-phase solitons become unstable in the whole first 
gap when Wq > 0.53. For comparison, out-of-phase soli- 
tons in the first gap do not exist in local FT symmet- 
ric defective superlattices with positive defects [2J|. It 
shows that the nonlocality expands the existing and sta- 
ble ranges of solitons. 

The field distributions of defect solitons are shown in 
Figs. [3l^a)-[3l^f) for different nonlocal degrees and prop- 
agation constants, which correspond to the cases rep- 
resented by circls in Fig. [UJa). Figures ^a.) and ^h) 
show that in-phase solitons exist in the semi-infinite gap, 
while solitons in the first gap are out-of-phase as shown 
in Figs. |3Ic)-12If). As the nonlocality degree increases, 
soliton power increases but the shape of their field dis- 
tribution changes very little. The propagations of soli- 
tons are also simulated based on Eqs. ([T]) and ([2]), and 
1% random- noise perturbations arc added into the initial 
input to verify the results of the linear stability analy- 
sis. The propagations corresponding to solitons in Figs. 
[3i;a)-l3i;f) are shown in Fig. |4i;a)-|4i;f). We can see that 
in-phase solitons are stable in the semi-infinite gap and 
out-of-phase solitons are stable in the first gap for the 
low propagation constants. 

Figures [S][7] show the results for zero defects (e = 0) . 
The properties of in-phase solitons in the semi-infinite 
gap and out-of-phase solitons in the first gap for zero 
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FIG. 5. (Color online) (a) Soliton Power P versus propagation 
constant b for different value of d. (b) Unstable growth rate 
Re{5) versus propagation constant b for out-of-phase solitons 
in the first gap. For all cases Vb = 6, Wo ~ 0.4, and e — 0.0 
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FIG. 6. (Color online) The complex fields (solid blue: 
real part; dotted red: imaginary part) and refractive-index 
changes (dashed black) for in-phase solitons in the semi- 
infinite gap at (a) b = 3.3, d = 1 and (b) b = 3.3, d — 3, 
and for out-of-phase solitons in the first gap at (c) b = 1.8, 
d = 1, (d) fe = 1.8, d = 3, (e) 6 = 2.2, d = 1 and (f) b = 2.2, 
d = 3, respectively. 



defects are almost the same as those for positive de- 
fects, except the positions of the cutoff points. From Fig. 
[SJa) we can see that the cutoff points approach the edge 
of the gap, so in-phase solitons are stable in the whole 
semi-infinite gap when Wq = 0.4. For comparison, the 
local defect solitons in the semi-infinite gap are unsta- 




FIG. 8. (Color online) (a) Soliton Power P vs propagation 
constant b for in-phase solitons in the semi-infinite gap and 
the first gap. (b)Cutoff point be versus the depth of defects, 
(c) Unstable growth rate Re{S) versus propagation constant b 
for the in-phase solitons in the semi-infinite gap. For all cases 
Vo = 6, W^o = 0.4, and e = -0.3 except for (b). 
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FIG. 7. (Color online) (a)-(f): Evolutions of defect solitons 
corresponding to Figs. [Hc)-[nii)) respectively. 
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FIG. 9. (Color online) (a)-(d) The complex fields (solid blue: 
real part; dotted red: imaginary part) and refractive-index 
changes (dashed black) for in-phase solitons in the semi- 
infinite gap at(a) 6 = 4, d = 0.5 and (b) 6 = 4, d = 1, or 
in the first gap at (c) b = 2.2, d = 0.5 and (d) b = 2.2, d = 1, 
respectively. 



ble near the edge of the gap[2J|. For nonlocal super- 
lattices, we find that in-phase solitons are always stable 
when Wo < 0.48. As increasing Wq, the stable region be- 
gins to shrink, and in-phase solitons are unstable in the 
whole semi-infinite gap when Wq > 0.76. 

Figure[5]Jb) is the perturbation growth rate Re{S) ver- 
sus propagation constant for out-of-phase solitons in the 
first gap with different nonlocality degree when Wq = 0.4. 
We can see that out-of-phase solitons are stable for the 
low propagation constant. As increasing Wq, the stability 
region of out-of-phase solitons vanishes after Wq > 0.51. 

The examples of in-phase solitons in the semi-infinite 
gap and out-of-phase solitons in the first gap for zero de- 
fects are shown in Figs. [ni[a)-[S]Jf), and the corresponding 
propagations are shown in Fig. [T] We can see that in- 
phase solitons in the semi-infinite gap have more peaks 
than those for positive defect. Figure [7] shows that in- 
phase solitons are stable at the edge of the gap whereas 
out-of-phase solitons are stable with low propagation con- 
stants. However, out-of-phase solitons in local FT sym- 
metric defective superlattices do not exist in the first gap 
for zero defects ! 
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FIG. 10. (Color online) (a)-(d): Evolutions of defect solitons 
corresponding to Figs. [Ha)-[8jd), respectively. 



For negative defects, we assume e = —0.3 and the re- 
sults are shown in Figs. [5]- [TUl We find that in-phase 
sohtons in the semi-infinite gap exist only for a weak 
nonlocality degree (d < 1.2), and our numerical program 
can not find any stationary solution for d > 1.2. The lin- 
ear stability analysis shows that in-phase solitons in the 
semi-infinite gap are unstable in their whole existence 
regions, as shown in Fig. ^c). In the first gap, differ- 
ing from positive and zero defects, in-phase solitons are 
found and stable in their whole existence regions when 
Wo < 0.56. As increasing Wq, in-phase solitons in the 
first gap become unstable for the large propagation con- 
stant when Wq < 0.78. Solitons are all unstable after 
Wo > 0.78. 

We can also see that there exists a cutoff point of 
propagation constant above which in-phase solitons exist. 
This feature is similar to the case for positive and zero 
defects. Figure [IJb) shows that the cutoff point shifts 
toward the lower propagation constant with increasing 
the depth of negative defects, and the existence region 
of in-phase solitons in the first gap increases too. When 
e = —0.435, the cutoff point arrives at the edge of the 
gap and then in-phase solitons exist stably in the whole 
first gap for Wq = 0.4. 

The examples of in-phase solitons in the semi-infinite 
gap and the first gap are shown in Figs. inUa)-|ntd), and 
the corresponding propagations are shown in Figs. llOf a)- 
[TUf d). One can see that in-phase solitons are unstable in 
the semi-infinite gap whereas they are stable in the first 
gap. 

Finally, we study the case of PT symmetric superlat- 
tices above the phase transition point. We find that in- 
phase solitons for positive and zero defects can exist sta- 
bly even if PT potentials are above the phase transition 
points when the value of Vq is small. We take Vq — 0.2 
and the phase transition point Woc = 0.05, as shown 
in the fig. [ijc). Figures [TlTa) and ITlT b) show the real 
and imaginary parts of bandgap structures for Vq = 0.2 
and Wq = 0.15, respectively. One can see that the cor- 
responding band diagram is partially complex. Figures. 
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FIG. 11. (a) Real and (b) imaginary parts of the bandgap 
structures for the superlattice with ei — 0.5, Vb ~ 0.2 and 
Wo = 0.15, respectively. 




FIG. 12. (Color online) The complex fields (solid blue: 
real part; dotted red: imaginary part) and refractive-index 
changes (dashed black) of solitons for (a) positive defect 
e = 0.3, (b) zero defect e = 0, and (c) negative defect 
e — —0.3. (d)-(f) Evolutions of defect solitons correspond- 
ing to these solitons in (a)-(c), respectively. For all cases 
Vo = 0.2, Wo = 0.15, b = 0.2, and d = 1. 

[T^ a)-fT^c) show the examples of solitons above the phase 
transition points for positive (e = 0.3), zero (e — 0), and 
negative defects (e = —0.3), respectively. Figures [T2l d)- 
n^ f) show the corresponding propagations of those soli- 
tons. Supported by the nonlocal nonlinearity (d = 1), de- 
fect solitons for positive and zero defects are stable when 
Wo is above the phase transition point, whereas they are 
unstable for negative defects. It is noted that the imag- 
inary parts of the complex fields for solitons shown in 
Fig. [T^ are still very small comparing with their real 
parts, although the PT potentials are above the phase 
transition points. 



IV. CONCLUSION 

We have studied the existence and stability of defect 
solitons supported by parity-time symmetric nonlocal su- 
perlattices. Unusual properties are found in this system. 
The nonlocality can drastically affect the existence and 
stability of defect solitons. For positive or zero defects, 
in-phase solitons can exist stably in the semi-infinite gap 
and out-of-phase solitons can exist stably in the first gap. 
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For negative defects, in-phase solitons are stable in the 
first gap, but unstable in the semi-infinite gap for the 
weak nonlocality degree. The values of Wq and Vq affect 
the soliton stability strongly. When Vq is large, defect 
solitons are unstable unless Wq is far below the phase 
transition point Woe- However, when Vo is small , the 
stable solitons can be found for positive and zero defects 
even if Wq is above the phase transition points. These 
properties of defect solitons in PT symmetric nonlocal 
superlattices are obviously different from those in the lo- 



cal PT superlattices. 
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